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> ■ 1 INTRODUCTION 
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j_j . In the highly successful model of dark matter-dominated hierar- 
C\3 chical galaxy formation, clusters of galaxies grow over time by 
accreting matter from their surroundings, with a well-defined dis- 
tribution of halo masses ranging from isolated gal a xies to large 



ABSTRACT 

We measure the diversity of galaxy groups and clusters with mass M > 10 13 /i~ 1 Af Q , in 
terms of the star formation history of their galaxy populations, for the purpose of constraining 
the mass scale at which environmentally-important processes play a role in galaxy evolution. 
We consider three different group catalogues, selected in different ways, with photometry and 
spectroscopy from the Sloan Digital Sky Survey. For each system we measure the fraction 
of passively-evolving galaxies within i?200 an d brighter than either M r = —18 (and with 
z < 0.05) or M r = —20 (and z < 0.1). We use the (u — g) and (r — i) galaxy colours to dis- 
tinguish between star-forming and passively-evolving galaxies. By considering the binomial 
distribution expected from the observed number of members in each cluster, we are able to 
either recover the intrinsic scatter in this fraction, or put robust 95% confidence upper-limits 
on its value. The intrinsic standard deviation in the fraction of passive galaxies is consistent 
with a small value of < 0.1 in most mass bins for all three samples. There is no strong trend 
with mass; even groups with M ~ 10 13 /i _1 Af Q are consistent with such a small, intrinsic dis- 
tribution. We compare these results with theoretical models of the accretion history to show 
that, if environment plays a role in transforming galaxies, such effects must occur first at mass 
scales far below that of rich clusters, at most M ~ 10 13 M Q . 
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tnbution ot nalo masses ranging irom isolated galaxies to large 
groups (e.g. lLacev & Cold 1 1994 Izhao et al.l 120031 ; iBerrier et all 
2009; iMcGee et alj|2009t) . If galaxy evolution is sensitive to the 
mass of the host dark matter halo, then differences in mass accre- 
tion history should be reflected in the population residing within 
groups and clusters JlvlcGee et alj2009h . For example, the increase 
in the fracti on of young galaxies in cluster cores as a funct ion of 
redshift (e.g. lButcher & Oemlell984lMargoniner et alj200lh can 
be linked to a corresponding increase in infall rate lEllingson et al.l 
l200ll ; lKodama & BoweJl200lUHaines et alj|2009l) . and the change 
in population with distance from the cluster centre can be used to 
infer a transformation timescale (e.g. Balogh et al. 20(3). 

Recently, IMcGee et al] fe009l) have shown that the cluster-to- 
cluster scatter in galaxy populations is a potentially powerful in- 
dicator of "pre-processing" — environmentally-driven transforma- 
tion that may have occurred in galaxies before they were accreted 
into t he current structure (e.g Zabludoff & Mulchaev 1998; Li et al] 
2009). For example, if the only environmental effect on galaxies oc- 
curs at mass scales well below that of rich clusters, one should ex- 



pect a very homogeneous cluster population. On the other hand, if 
galaxies are transformed only once they are accreted into very mas- 
sive systems, than the stochasticity of this process will lead to more 
diversity in present-day cluster populations. In fact, s uch an argu- 
ment was recently invoked by IPoggianti et al] ]2006h . who claim 
that an increase in diversity in groups with velocity dispersions be- 
low a ~ 500km/s is evidence that transformations are likely oc- 
curring at those scales. However, the uncertainties on the fraction 
of galaxies with [Oil] emission in a single system are large, and 
they increase with decreasing mass. Thus it is not clear whether or 
not the intrinsic scatter is such a strong function of mass. 

With the advent of very large surveys of nearby galaxies, the 
populations of galaxy groups and clusters and their correlations 
with mass, X-ray p r operti e s and redshift are fairly w ell established 
(e.g. Fairlev et al. 2002; De Propris et al. 2004; Popesso et al 
I2004L IWake et al] |20 05; Wei nmann etafl 120061; [Aguerri et al 



| 2007l: iFinn et al] 120081; iKimm et al] 120091 ; lHansen et al] 120091 : 
iBarkhouse et alj2009l ; lLu et al.l2009T) . However, relatively little at- 
tention has been giv en to the varia t ion bet ween clusters. One im- 
portant exception is Popes so~et al] d2007bh . who analyzed 79 X- 
ray clusters in the Sloan Digital Sky Survey (SDSS), and find a 
r.m.s. scatter of 0.19 in the fraction of blue galaxies. However, they 
do not have a large enough sample to study this as a function of 
cluster mass, which is crucial for identifying a putative transfer- 
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mation scale. Furthermore, the measured r.m.s. includes a contribu- 
tion from the statistical uncertainties on individual measurements, 
and thus is an uppe r limit on the intrinsic variation. Very recently, 
lHaines et alj J2009h use Spitzer data of 30 X-ray luminous clusters 
at z < 0.3 to measure the fraction of strongly star-forming galax- 
ies (> 8M0yr _1 ), and find a remarkably small intrinsic scatter, 
consistent with zero once the trend with redshift is accounted for. 

In this paper, we revisit the issue using three different samples 
of galaxy clus ters, based on the SPS S. O ur results are consisten t 
with those of IPoggianti et all d2006h and IPopesso et alj b007bh . 
but we take the extra step of measuring the intrinsic cluster-to- 
cluster variation, an d comparing this with the model predictions of 
iMcGee et alj d2009h - Throughout the paper we adopt a cosmology 
with Q m = 0.3, Aa = 0.7, and parameterize the Hubble constant 
as Ho = 100/ikm s _1 Mpc _1 . 



2 SAMPLE SELECTION 

For our purposes, we require a large, homogeneous sample of 
galaxy groups and clusters, together with a simple, reliable and 
sensiti ve measurement of the galaxy population within them. The 
SDSS jYorketal.ll2000h , with its large size, homogeneous data, 
and highly complete spectroscopy, is particularly well-suited to this 
type of study. We use data from the DR6 spectrosco pic sample, 
making use of the NYU-VAGC of iBlanton etail J2005h . The data 
are unbiased for r < 17.77, and colours and luminosities are mea- 
sured from Petrosian magnitudes, k- corrected to z = 0.1 using 
KCORRECT telanton & Row eis 2007). We make no correction for 
spectroscopic completeness; this is a s mall correction with l ittle 
dependence on luminosity or colour (e.g. lPopesso et alj2007bl) . so 
has no impact on these results. 

There are many different ways to find galaxy clusters, usu- 
ally based on either galaxy position (with or without redshifts, and 
possibly using colour information) or X-ray emission from the in- 
tracluster plasma. The ideal sample is complete, uncontaminated, 
and has a reliable, observable property that can be related to the 
total mass of the system. For our purposes, completeness is per- 
haps the most critical, as this allows us to put robust upper limits 
on the scatter we observe. We will use three different cluster cat- 
alogues for this analy sis. The first, wh ich we call the "Halo" cata- 
logue, is taken from lYang et alj l |2005t) , who use some prior infor- 
mation from theory and observation to associate every galaxy with 
a "halo". The mass is then given by rank-ordering the groups by 
their total luminosity or stellar mass (we use the latter), and associ- 
ating that with a theoretical dark matter mass function. The secon d 
catalogue, also based on optical data, is from lBerlind et alj {2006). 
We refer to this as the "FOF" catalogue, as it is a traditional friends- 
of-friends algorithm, with parameters calibrated to match numeri- 
cal simulations. Masses for these systems are calculated from their 
velocity dispersions, assuming virial equilibrium. Finally, we con- 
sider an X-r ay selected sample of cluste rs, which consist of all HI- 
FLUGCS teeiprich & Bohringenl2002l) clusters in SDSS, and C4 
(Miller etalj|2005h clusters cross-correlated with the ROSAT cata- 
logue. For cl usters in the HIFLUGCS sam ple, we use the masses 
measured by Reiprich & Bohringer ( 2002), under the assumption 
of hydrostatic, isothermal gas; these have typical uncertainties of 
1 0-30 per cent. For the C4 clu sters, we use the M — L x relation 
of iReiprich & Bohringej d2002h to convert ROSAT luminosities to 
mass. 

Galaxy properties are known to correlate strongly with both 
clustercentric radius, and the limiting luminosity or stellar mass 



of the galaxy sample Pope sso et alj (e.g. l2007bl) . We will select 
all galaxiesjj within the estimated projected virial radius r2oo, 



determined from the cluster mas£J. This is done in different 
ways for each sample . The X -ray sample uses r2oo computed by 
IReiprich & Bohringej j2002h fo r the HIFLUGCS clusters, and the 
results of IPopesso et alj (2007a), based on fitting a King profile to 
the galaxy distribution, for the RO SAT-matched C4 clusters. For 
the FOF group catalogue we calculate r2oo from the velocity dis- 
persions, using eq uation 8 in | Finn et alj J2005I) . The masses of the 
Halo clusters from lYang et alj J2007h are measured by ranking the 
clusters by stellar mass and comparing to dark matter mass func- 
tions, but the mass calculated is for a radius with a matter overden- 
sity of 1 80 which, for their choice of cosmology, corresponds to 43 
times the critical density. Therefore their masses and radii are about 
2.2 times larger than A/200 and -R200, assuming a flat potential. We 
apply this correction to make them more comparable to the other 
two samples. However, we need not concern ourselves overmuch 
with ensuring complete consistency between samples; indeed, dif- 
ferences between them help show explicitly that our main conclu- 
sions are completely independent of these details. Both statistical 
and systematic u ncertainties (e.g. proje ction of field galaxies along 
the line of sight. IPopesso et alj|2007bh associated with this mea- 
surement will only lead to increase the scatter we measure; thus 
our upper limits will still be robust. 

We will select galaxies based on luminosity, and will consider 
two versions of each sample. The main catalogue is limited to z < 
0.05, and we use all galaxies brighter than M r + 5 log h — — 18; 
for comparison we consider a "bright" sample consisting of M r + 
5 log h < —20 galaxies, and including clusters out to z = 0.1. Our 
results on the scatter are comparable for both samples. 



3 DEFINITION OF PASSIVE GALAXIES 

We will take advantage of the well known fact that galaxies appear 
to primarily divide into two classes. One has predominantly red 
colours, and little star formation , while the other consists of blue 
and actively forming sta r s (e.g. Istrateva et al]|200ll ; iBaldrv et alj 
120041 ; iKimm et alj 120091 : IPozzetti et alj I2009L and many others). 
Most simply, this distinction can be made using a colour-magnitude 
diagram to isolat e red-sequence and blue-cloud galaxies (e.g. 
IBaldrv et alj|2004l) . However it is also known that a single optical 
colour cannot distinguish betwee n dusty, star-for ming galaxies, and 
truly old, passive gala xies (e.g. IPopesso et aul2007bo . However, 
Wolfet al. (2005, 2009) have shown that this distinction can be reli- 
ably made using two colours: one bracketing the 4000 A break, and 
another at longer wavelengths. We show this for our full sample of 
SDSS galaxies, in Figure [T] each panel shows the (it — g) colour 
as a function of (r — i), for galaxies within 0.1 magnitude of the 
r— band luminosity shown (all colours are k-corrected to z = 0.1). 
The existence of two populations is remarkably clear, as is the fact 
that a cut in (u — g) colour alone would include galaxies from both 
populations. 

We choose to select passive galaxies as those within the ellipse 
shown in Figure Q] The centre of the ellipse is a smooth function 
of Mr, but the orientation and size is kept fixed. These ellipses 



1 Only clusters with full SDSS coverage at r < R200 are included in our 
analysis. 

2 We define -R200 t0 be the radius within which the total matter density is 
200 times the critical density of the Universe. 
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Figure 1. Colour-colour diagrams (k-corrected to z=0. 1) for a random sub- 
set of SDSS galaxies in our sample within 0. 1 mag of the M r magnitude 
indicated in the top-left corner (h = 1). The galaxies are clearly separated 
into two populations, and we associate those galaxies within the red ellipses 
as "passive" galaxies, with little or no ongoing star formation. The ellipse 
is fixed in size and orientation, but the centre is a linear function of M r . 

are not optimized in any rigorous way, but reasonable variations in 
definition have no influence on our results. 

For comparison we will also measure a simple "red fraction", 
based on the (it — r) colour as a function of M r . We fit the red 
sequence with a slope of —0.0818, and a colour (u — r) — 2.8 at 
M T + 5 log h — —20, and choose red galaxies to be those that are 
up to 0.25 mag bluer than this line. Our results on the scatter of the 
red fraction are completely insensitive to the details of this choice. 



4 RESULTS 

In Figure[2]we show the fraction of passive galaxies in every clus- 
ter in the three catalogues, as a function of its mass, for the "main" 
sample (M r + 51og/i < —18). Uncertainties on the fractio ns are 
computed using the full binomial distribution (Gehrels 198^). We 
see that both the Halo and X-ray samples of clusters are in good 
agreement where they overlap, despite the different selection crite- 
ria, and definitions of mass and radius. The passive fractions are 
systematically larger than the M r + 5\ogh < —18 sample as 
a whole (31%, shown as the horizontal line), on average. There 
is a weak trend with mass, such that the most massive clusters 
have larger passive fractions, on average. The passive fractions 
are systematically lower than the red fractions based only on the 
(u — r) colour (top panel), due to the contamination from dusty- 
spiral galaxies in the latter. However, the scatter from cluster-to- 
cluster is similar whether we consider the red or truly passive frac- 
tion. 

The scatter in Figure [2] becomes very large at lower masses. 
This is very similar to the trends shown in Figs. 4b and Fig.6 of 
Pogaianti et al. (2006), where they find the scatter in the fraction 
of galaxies with [Oil] -emission increases sharply in systems with 
velocity dispersion below a ~ 500 km/s (M ~lx 10 14 /i _1 M Q ). 
However, the uncertainties on individual measurements are also 



T 




aasffico^ □ □ 



1QU IQU [QUS 

Cluster mass (h _1 M G ) 

Figure 2. Bottom panel: The fracti on of passive gal axies in each clus- 
ter within the o ptical Halo catalogu e lYang et al]|2005l open squares), the 
FOF catalogue i Berlind et al. 20061 open triangles) and the X-ray matched 
catalogue (filled circles). In all cases we use the main catalogue, with 
M r + 51ogh < —18. Passive galaxies are defined as those within the 
elliptical areas in the colour-colour diagram of Figure [T] The dotted line 
shows the passive fraction in the full SDSS sample (31%), considering 
M r + 5 log h < — 18 galaxies. Error bars are 68% confidence limits using 
the binomial probability distribution, and are shown only on the X-ray sam- 
ple, for clarity; the uncertainties on individual points at low masses get very 
large. Top panel: Similar, but here we show the fraction of red galaxies, 
defined using only the (u — r) colour. 

much larger in these lower mass systems, since they have fewer 
members. Our task now is to recover the intrinsic scatter from these 
observations. 

In Figure[3]we show the average passive fraction as a function 
of mass, for both samples. The black "error bars" show the standard 
deviation in this fraction (they are not the error on the mean, which 
is much smaller). The red "error bars" show the standard deviation 
expected from statistical uncertainties alone, assuming no intrinsic 
variation. To compute this we treat the observed fraction for each 
cluster as randomly drawn from a binomial probability distribution 
defined by the number of galaxies in the cluster and with an expec- 
tation value given by the mean value in each mass bin. We neglect 
any error in the mean value itself, which is small. 

This demonstrates that the variance in passive fraction is not 
much larger than expected due to the statistical uncertainties alone. 
The close agreement between all three samples also suggests that 
systematic uncertainties related to sample selection, and mass or ra- 
dius definition, contribute much less to the observed variance than 
the statistical uncertainties. 

For each mass bin we simulate the binomial distribution 10000 
times, and compute the variance of each simulation. By adding in 
quadrature an assumed amount of intrinsic (Gaussian) variation, 
ranging from to 0.25, we estimate the 5% and 95% confidence 
limits on this intrinsic scatter, and specifically the probability that 
it is non-zero. When the latter is greater than 95%, we measure the 
instrinsic scatter as the quadrature difference between the observed 
and predicted statistical standard deviations; otherwise, we report 
an upper limit. In Figure|4]we show these measurements and upper 
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Figure 3. Top panel: The points show the average passive fraction among 
all clusters, in bins of equal mass; only points with at least five contributing 
clusters are shown. The filled circles repre sent the X-ray cat alogue, while 
the open squares are the Halo cata logue o f Yang et al. (2005) and the open 
triangles are the FOF catalogue of lBerlind et alj (2006). Black ' 'error bars" 
indicate the standard deviation within each bin (not the error on the mean). 
Red bars represent the standard deviation expected due to binomial sam- 
pling fluctuations alone, assuming no intrinsic scatter. The horizontal, dot- 
ted line indicates the passive fraction in the whole Mr < — 18 SDSS sam- 
ple, 31%. The observed variance is typically very similar to that expected 
from the statistical uncertainties alone. Bottom panel: The same, for the 
M r + 5 log h < — 20 sample, which has more clusters but fewer galaxies 
per cluster. Here the horizontal, dotted line indicates the passive fraction in 
the brighter M r < -20 SDSS sample, 49%. 



limits on the scatter, as a function of mass. Generally, the scatter is 
remarkably small, consistent with 0.1 or less. Importantly there is 
at most a weak trend with cluster mass: even at M ~ 10 13 /i _1 Mq 
the intrinsic variation is small. The fact that a large scatter is ob- 
served in these systems is due to the fact that the uncertainties as- 
sociated with an individual group are large, limited by the number 
of members. With a large enough sample, however, it is possible to 
overcome these uncertainties and allows us to put interesting limits 
on the intrinsic scatter. 



5 DISCUSSION AND CONCLUSIONS 

We have shown that the fraction of passive galaxies in clusters has 
remarkably small scatter between systems, generally 0.1 or less. 
We ca n compare this directly with the predictio ns of iMcGee et al.l 
(2009), who use the galaxy formation model of [Font et al . (2008) 
to calculate the rate at which galaxies are accreted into groups 
and clusters, and the amount of time they spend wit hin haloes of 
a give n mass during their history. More specifically, McG ee et al.l 
( 2009) show the fraction of galaxies that have been inside a halo of 
mass M > M trU nc for a time t > Ttrunc- Under the assumption 
that such galaxies are the only ones that turn passive, we can com- 
pare those predictions with our results. Although these predictions 
make use of a semi-analytic model to trace galaxies, the results are 
determined primarily by the dark matter growth history, which is 
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Figure 4. We show the intrinsic la r.m.s scatter on the passive fraction 
from cluster to cluster, in mass bins with at least five clusters. The X-ray 
clusters are indicated with filled circles; the Halo clusters with open squares, 
and the FOF clusters with open triangles. Where the observed variance is 
significantly larger than the statistical uncertainty, we have represented the 
95% confidence limits with error bars. In cases where we are unable to 
measure a significant intrinsic scatter, in light of the statistical uncertainties, 
we have placed the point at the 95% confid ence upper limit, in dicated with 
a downward arrow. Model predictions from McGee et al. ( 2009) are shown 
as the lines. Red, green and blue lines indicate log(hMtrunc) = 14, 13 
and 1 2, respectively. In each case the solid line represents Ttrunc = 3Gyr, 
with the dotted lines representing 1 Gyr (lower scatter) and 6 Gyr (higher 
scatter). Top panel: The main cluster samples, used throughout most of this 
sample, restricted to M r + 51og/i < —18 and z < 0.05. Bottom panel: 
The larger, "bright" cluster sample, restricted to M r + 5 log h < — 20 and 
z < 0.1. 



taken from the Millennium simulation jSpringel et al. 1 120051 ). and 
are largely undependent of the recipes used to mode l the galaxies. 

We plot predictions from IMcGee et alj J2009h together with 
our data, in Figure[4] Note that the observed scatter will be affected 
by other systematic effects, such as field contamination and un- 
certainties in the mass and radius; thus all of our points are best 
thought of us strict upper limits on the true cluster population, 
which is what is being modeled. Lines are shown for three differ- 
ent Aitrunc and three different Ttrunc, as described in the caption. 
It is interesting first of all that the scatter predicted by the mod- 
els is quite comparable to that observed, even in the lowest mass 
groups; there is no evidence that the variance in group properties 
significantly exceeds what would be expected from these simple 
infall-based models. The fact that these groups have passive frac- 
tions significantly greater than the global average, with a small in- 
trinsic scatter, implies that A/trunc < 10 13 /i _1 Mq. That is, star 
formation in galaxies must be shutting down long before they enter 
cluster-sized haloes. Treating all the points as strict upper limits, 
even the scatter in the most massive systems puts interesting limits 
on these parameters. If star formation is only truncated upon accre- 
tion into clusters, A/trunc > lO 14 /i -1 Af0, then the timescale must 
be T < lGyr to be consistent with all the cluster samples presented 
here. Such a model, however, would predict greater scatter than ob- 
served at M — 2 x 10 14 h" 1 Mq, and moreover that systems with 
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M < Mtrunc should resemble the field population, which is ruled 
out. 

The results of lHaines et ID 12009b potentially provide even 
stronger constrai nts. Their sample of clusters with ma sses of ~ 2- 
20 x 1O 14 M dOkabe et alj|2009t IZhang et alj |2008) has an rms 
of 0.03 after accounting for a trend with redshift, and is consistent 
with no intrinsic scatter. At face value, this suggests that not only 
must Aftrunc be low, but T trun c must also be small, < 1 Gyr. How- 
ever, their star formation rate threshold is quite high, and this is 
probably why the average fraction of passive galaxies at the low- 
est redshifts in their sample is 0.95, substantially higher than in the 
present sample. This, together with the wider redshift range cov- 
ered, makes it difficult to directly compare with our results. 

In conclusion, the fraction of passive galaxies and the variance 
in this fraction from system to system suggests that star formation is 
shut off in galaxies within groups with masses 
thus, "pre-processing" is crucial to explain the observed properties 
of today's clusters. These constraints can be significantly improved 
by increasing survey depth, so there are more members per clus- 
ter, and by increasing the volume surveyed so there are more con- 
tributing clusters. Mass measurements with smaller statistical un- 
certainty (for example from X-ray observations with resolved tem- 
perature profiles) would also be helpful in reducing the scatter, as 
there is a small trend for the passive fraction to increase with mass. 
Perhaps most importantly, repeating the analysis at higher redshift 
will be valuable, as the scatter associated with large Ttrunc mod- 
els is predicted to increase significantly with redshift (McG ee~et al.l 
l2009h . 
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